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We use a path integral quantum Monte Carlo method to simulate excitons and biexcitons in core 
shell nanocrystals with Type-I, II and quasi-Type II band alignments. Quantum Monte Carlo tech- 
niques allow for all quantum correlations to be included when determining the thermal ground state, 
thus producing accurate predictions of biexciton binding. These subtle quantum correlations are 
found to cause the biexciton to be binding with Type-I carrier localization and strongly anti-binding 
with Type-II carrier localization, in agreement with experiment for both core shell nanocrystals and 
dot in rod nanocrystal structures. Simple treatments based on perturbative approaches are shown 
to miss this important transition in the biexciton binding. Understanding these correlations offers 
prospects to engineer strong biexciton anti-binding which is crucial to the design of nanocrystals for 
single exciton lasing applications. 



I. INTRODUCTION 

The properties of excitons in semiconductor quantum 
dots have in recent years been the focus of consider- 
able study, both as a source of new fundamental physics 
and as a quantum electronic system that can be en- 
gineered for a range of applications including optical 
sources and detectors. A striking example of such en- 
gineered quantum confinement is found in Type-II core- 
shell nanocrystals^^ as depicted in Fig. [I] Exploiting in- 
dependent control of the electron and hole wavefunctions 
it is possible to tune the exciton energies, inter-particle 
Coulomb and exchange interactions, and transition life- 
times to a far greater degree than in Type-I band-aligned 
structures. Of particular practical significance is the abil- 
ity to generate positive exciton-exciton (X-X) interaction 
energies (anti-binding) as a result of such spatial manip- 
ulation of electrons and holes. Then the photon energy 
required to generate a biexciton in a quantum dot al- 
ready containing an exciton is greater than the exciton 
recombination energy. This is an important pre-requisitc 
for achieving lasing in the single exciton regime since an 
incident photon with energy resonant with the exciton 
energy may stimulate emission but cannot be absorbed. 
Exciton-exciton interaction energies of up to +110 meV 
have been reported in CdS/ZnSe core-shell structures^ 
but single exciton lasing from nanocrystals has yet to be 
achieved. Recently 'dot-in-rod' structures in which an 
approximately spherical CdSe nanocrystal is embedded 
inside a CdS rod have also been shown to possess control- 
lable quasi- Type II behaviorJSHni anc j thus are candidate 
systems for single exciton lasing. 

Previous calculations of the exciton-exciton interaction 
energy (A^ = E xx — 2E X where E x and E xx are the ex- 
citon and biexciton total energies) in Type- II nanocrys- 
tals have mostly relied on first or der pertu rbation theory 
and assumed spherical symmetry! 3 * 4 * 12 * 1 ^ These calcula- 
tions implicitly assume the limit of strong confinement, in 
which the energetic separation of the single particle states 
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FIG. 1. CdTe/CdSe Type-II core/shell nano-crystal 
schematic and band edges, with electron and hole probabil- 
ity densities within a biexciton. a) 1.95 nm core radius and 
0.25 nm shell thickness, b) 1.95 nm core radius and 2.5 nm 
shell thickness. 



is much greater than the inter-particle interaction energy. 
They provide a first approximation to the increase in 
carrier repulsion upon growth of a Type-II aligned shell 
layer, but consistently overestimate the repulsive effect 
since they do not include any spatial correlations in the 
carrier wavefunctions. Measurements of A xx show that 
it is always negative (binding) for core-only nanocrys- 
tals, and becomes positive (anti-binding) only as the shell 
thickness is increased beyond a threshold valueP Per- 
turbative models miss this important transition between 
binding and anti-binding as the nanocrystal confinement 
passes from the Type-I to the Type-II regime. Korkusin- 
ski et at employed a configuration interaction (CI) ap- 
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proach using a tight binding basis set to calculate the ex- 
citon and biexciton binding energies in wurtzite core-only 
CdSe nanocrystals,™ finding that biexciton anti-binding 
occurs for nanocrystals smaller than about 4 nm since the 
crystal field provides greater localization of the quantum 
confined hole than the electron. In many situations CI 
calculations converge only slowly as the basis size is in- 
creased. This is because the influence of higher lying 
states drops off as 1 j AE where AE is the energy of the 
basis state relative to the system energy, and the den- 
sity of states tends to increase rapidly as the radius in- 
creases and as higher energy states are included. These 
two trends combined make it computationally challeng- 
ing to include sufficient states to be confident of good 
convergence. 

In this work we report on path integral quantum Monte 
Carlo (PI-QMC) calculations of exciton binding ener- 
gies and exciton-exciton interaction energies, with a full 
treatment of the quantum correlations in the many body 
system. PI-QMC does not suffer from the convergence 
problems of CI calculations mentioned above. The QMC 
approach has been used previously to calculate the prop- 
erties of simple nanocrystals,™ the results of which illus- 
trated the importance of correlation effects in determin- 
ing the biexciton binding, but have yet to be applied to 
anti-binding scenarios and to multishell heterostructures. 
Unlike pcrturbative approaches, our results capture the 
experimentally observed transition from binding to anti- 
binding character in Type-II nanocrystals. We use this to 
study the Type-II regime in CdTe/CdSe and CdS/ZnSe 
core/shell nanocrystals, and make comparisons to pre- 
vious experiemental results. We then invetigate an in- 
verted Type-I core/shell structure, ZnSe/CdSe, and a 
quasi- Type II dot-in-rod structure, CdSe/CdS, both of 
which have been suggested as possible candidates for 
strong anti-binding. 



II. MODEL 

We focus on the binding/anti-binding transition of the 
biexciton where an accurate treatment of the correlation 
energy is crucial to obtain the correct magnitude and 
sign of the biexciton binding.™ As discussed in Ref. IT51 
we choose a simplified single-band effective mass model 
to enable essentially exact determination of correlation 
energy in the binding transition, while losing some of the 
details of atomistic models. We model the nanostructures 
using a biexciton Hamiltonian of the form, 
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where the kinetic energy arises from parabolic bands, 



and there are separate electron and hole confining poten- 
tials, V e and V/j, arising from the band edges in the core 
shell nanocrystal. 

V dot = V e (r ei ) + V e (r e2 ) + V h {r hl ) + V h (r h2 ). (4) 

The exciton Hamiltonian is a simple reduction from this 
form. This Hamiltonian treats carrier propagation in the 
nanocrystal systems within the single-band effective mass 
approximation, and we model the heterointerfaces as 
step-like potentials in the conduction and valence bands. 
This simplification of the semiconductor band structure 
is the main limitation of the present model, and one 
might expect more accurate predictions to be achieved 
using a multi-band description. The latter would be ex- 
pected to increase somewhat the degree of correlation 
due to Coulomb interactions, as it is known to reduce the 
energy spacing between quantum confined valence band 
states as a result of mixing effects™ 

Our model uses a finite potential barrier for the sur- 
rounding matrix, and we assume a uniform dielectric con- 
stant throughout. We do not consider dielectric polariza- 
tion effects in this work.™ 1 To check that this approach 
does not introduce substantial errors, we performed a se- 
ries of calculations that included the dielectric self-energy 
of the carriers, by adding the self-energy potential to 
the confinement potential before running the PI-QMC 
algorithm. The self-energy potential as a function of ra- 
dial coordinate was calculated according to the work in 
Refs. [16] and [TTJ Our results confirmed that although 
inclusion of this potential affected the single particle en- 
ergy levels significantly, it left the exciton binding and 
exciton-exciton interaction energies unchanged, in agree- 
ment with previous work™ 

We consider several different types of nanocrystals, ini- 
tially Type-II core/shell structures for both electron/hole 
and hole/electron confinement. Later we address some 
suggestions about the possibility of large X-X interac- 
tions in inverted Type-I structures. Finally we present 
some results for core/rod nanocrystal structures, in 
which there is currently significant interest due to their 
excellent as-grown uniformity and very high quantum 
yields P 
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III. METHODS 



The interaction potential includes all pair-wise Coulomb 
interactions, using a uniform dielectric constant (dis- 



The PI-QMC method allows essentially exact treat- 
ment of correlation energy, with no basis-set or varia- 
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tional bias. It is based around a stochastic sampling of 
the many body thermal density matrix, 

(0 (R,R';^) = i(R|e-^|R / ), (5) 

where /3 = l/fe^T is the inverse temperature, R = 
(r ei , . . . , Th 2 ) are the particle coordinates, and the par- 
tition function, Z = ti(e~^ H ), normalizes the density 
matrix. 

Averages of any physical observables O can then be 
calculated by, 



(6)= / dRdR>(R,R';/3)(R|d|R') 
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In the imaginary time path integral method, the ther- 
mal density matrix is expanded, using the primitive 
approximation,^ into N slices, 
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with m the mass of the path's particle, and with a time 
step At = f3h/N. The 1/r singularity of the Coulomb 
interaction requires extra care, so we replace the Ar V^oui 
action term with the pair approximation to the Coulomb 
action.^ We then use the well known Metropolis Monte 
Carlo algorithm for the stochastic sampling of the dis- 
cretized thermal density matrix, in which a closed quan- 
tum path (such that Ro = Rjv) representative of a par- 
ticle is randomly walked through real space. 

The algorithm is run until the required accuracy is 
reached, which must be high, as the binding energy is 
the difference between large total energies. Thus these 
total energies must have a small absolute errors (typically 
±0.5 meV, unless otherwise indicated with error bars) so 
as not to result in binding energies with large relative er- 
rors. The PI-QMC method has several advantages: the 
use of the thermal density matrix naturally provides fi- 
nite temperature simulations, where the temperature is a 
controllable parameter; it requires no basis set informa- 
tion or trial data; and it can treat spatially complicated 
potentials.^ Finally, and most importantly for the work 
presented here, it treats the many body Coulomb corre- 
lations exactly. 



IV. TYPE II CdTe/CdSe 

We first study Type-II CdTe/CdSe quantum dots as 
described by Oron et alW who have presented experi- 
mental results for a such core/shell structures, in which 
a transition from the binding to anti-binding regime in 
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FIG. 2. Exciton-exciton interaction energy, A xx , 
shell thickness for a CdTe core diameter of 3.9 
mental data from Oron et a/P are black crosses 
imental uncertainty approximately given by the 
Blue diamonds show the PI-QMC results, and 
guide to the eye. The dashed blue line shows 
theory results. Inset shows the radial form of 
ment potential. 
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the X-X interaction energy is clearly visible as the shell 
thickness is increased. In our model of this CdTe/CdSe 
(core/shell) quantum dot, the potential minimum for 
holes lies in the core and for electrons lies in the shell. 
The band gaps for CdTe and CdSe are taken as 1.475 eV 
and 1.75 eV wit h the CdTe/CdSe valance band offset 
taken as -0.57 eV! 20 * 21 l The external potential band gap 
is assumed to be 4.832 eV with a valence band offset of 
1.325 eV. We take the electron and hole effective masses 
isotropically as 0.13m e and 0.35 m e respectively!^! All 
our simulations for this dot were carried out at 300 K. A 
dielectric constant of 6.65 is used throughout the simu- 
lation unit cell. 

Fig. [2] shows the exciton-exciton interaction energy, 
A xx , plotted against shell thickness. Experimental data 
points are taken from Oron et al)^ for a 3.9 nm diame- 
ter core. We see excellent quantitative agreement with 
these experimental results, demonstrating the impor- 
tance of a correct treatment of correlation in calculat- 
ing the exciton-exciton interaction energies, even in this 
rather strongly confined nanostructure. At the smallest 
shell thicknesses in Fig. [2] the modeled shell thickness is 
smaller than the lattice spacing and, as in the experi- 
ments, reflects an ensemble average over many nanocrys- 
tals. 

In Fig. [3] we plot A xx against shell thickness for a va- 
riety of core diameters. The same trend is evident for all 
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FIG. 3. Exciton-exciton interaction energy (A xx ) versus CdSe 
shell thickness for various different CdTe core diameters, as 
indicated. Lines are a guide to the eye. Inset shows the radial 
form of the confinement potential. 



FIG. 4. Exciton interaction energies, A^, as a function of 
CdSe shell thickness for various CdTe core diameters, as in- 
dicated. Perturbation theory result for 3.9 nm core are shown 
as a dashed blue line. Lines are a guide to the eye. Inset 
shows the radial form of the confinement potential. 



four data sets; for small shell thicknesses the biexciton is 
strongly bound, with large binding energies. Small shell 
thicknesses result in a quasi Type-I structure, where the 
shell is not thick enough to induce localization of the elec- 
tron, thus the electron is spread across the core, close to 
the hole, leading to binding of the biexciton. For thicker 
shells, localization of the electron in the shell begins, and 
the transition to positive A xx is seen, due to the increased 
electron-hole seperation and hence a reduced attractive 
Coulombic interaction energy, whilst the large repulsive 
hole-hole interaction from the core confined holes remains 
relatively unchanged. 

Core/shell nanocrystals with smaller cores have a 
smaller binding energy and a larger anti-binding energy, 
due to the stronger confinement of the holes in the core, 
giving a larger repulsive element to the binding ener- 
gies. The smaller the core size, the larger the final anti- 
binding, as the separation between the holes within the 
core remains small despite their mutual repulsion. 

We now compare our PI-QMC results to that of first 
order perturbation theory to assess the role of correla- 
tion. Quite good agreement is observed between the 
exciton interaction energies, A x , calculated by the two 
methods as shown in Fig. [3] A x is the energy change at- 
tributable to the Coulomb interaction, i.e., the negative 
of the binding energy. The Coulomb correlation between 
the electron and hole wavefunctions that are included 
in the PI-QMC results increase the binding energy by 
7-10 meV compared with the perturbation theory predic- 
tions. For the biexciton (Fig. [2| the difference in A xx 
between the two approaches is more marked. The per- 



turbative approach fails to predict any transition between 
negative and positive values of A^, and further, over es- 
timates the anti-binding significantly. The correlations 
included in the PI-QMC approach lead to a reduction 
in A xx of about 17meV in the core-only Type-I struc- 
ture, and which gradually increases to about 29 meV for 
the fully Type-II structure with a 2.5 nm shell. The in- 
creasing importance of the correlations with increasing 
shell thickness can be attributed to the departure from 
the strong confinement limit as the electron becomes less 
localized. Calculations performed for a 3.9 nm core with 
thicker shells (not shown) reveal that A xx saturates to 
about 70 meV at a shell thickness of around 7.5 nm. 

The role of correlations in the biexcitons can be exam- 
ined directly through the conditional probability density 
function for the electrons and holes as shown in Fig. [5] 
Each sub-panel A1-A4 in Fig.[5]shows a slice of the condi- 
tional probability density through the x-y plane of a 6 nm 
CdTe core diameter colloidal dot with a 2.5 nm CdSe 
thick shell. This geometry places it strongly in the anti- 
binding regime, as seen in Fig. [3j The small rectangle 
in each panel indicates the location of a small volume, 
1Z, which subtends a solid angle of 0.5 degrees between 
the two radii at either end of the rectangle. (The rect- 
angles are drawn wider than 0.5 degrees for illustrative 
purposes.) 

The conditional probability density in each panel along 
the row (l)-4)) is defined by 
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FIG. 5. Conditional probability densities are shown for a 6nm CdTe core diameter and 2.5 nm CdSe shell thickness in row 
A, for a 6nm CdTe core diameter and 0.25 nm CdSe shell thickness in row B and for a 2nm CdTe core diameter and 2.5 nm 
CdSe shell thickness in row C. The radial form of the confinement potential for each is illustrated. Shown in column 1) is 
g ee , a conditional electron (falling within the blue rectangle) and the resulting electron distribution, column 2) shows g e h, a 
conditional electron and resulting hole distribution, column 3) shows gu e , a conditional hole (falling within the red rectangle) 
and resulting electron distribution and column 4) shows guh, a conditional hole and resulting hole distribution. 



where i denotes the particle that must fall into the region 
1Z in order for the location of particle j to be sampled. 
Hence plotted in sub-panel Al is the probability of ob- 
serving an electron at a given position if a first electron is 
found in the region 1Z. Similarly Fig.[5]A2 shows the pair 
correlation density for holes when an electron is present 
inside the rectangle. Fig. [5] A3 and A4 are analogous 
plots for when a hole is located in the region U. In Fig. |] 
Al it can be seen that the electrons repel one another 
and sit on opposite sides of the dot; the same behavior is 
seen between the two holes in A4 with the holes sitting 
on opposite sides of the core. In A2 and A3 opposite 
charges can be seen correlated and being attracted to- 
wards the opposite charge located in the region 1Z. It 
is clear that in addition to the radial correlations, angu- 
lar correlations also play a role in reducing the electron- 
electron and hole-hole interactions to further reduce the 
total Coulomb energy in the system. 

We contrast this with the case of the bound biexciton, 
for a 6nm CdTe core diameter and 0.25 nm CdSe shell 



thickness shown in the second row of Fig. [5] again similar 
electron-electron and hole-hole repulsion is seen in Bl 
and B4. However, the strong quantum confinement of 
the electrons and holes keeps them well confined to the 
core, leading to the strong electron-hole interactions seen 
in B2 and B3. This results in a much increased overlap 
between the electron and holes charge densities and this 
increase in the attractive electron-hole interaction gives 
rise to a bound biexciton. The correlations indicate a 
state in which the electron and hole motions are strongly 
overlapping. 

In a dot with a small 2 nm diameter CdTe core and 
2.5 nm CdSe shell thickness, we see a significantly in- 
creased anti-binding. In the final row of Fig. [5] C2 and 
C4 shows the hole being strongly confined in the small 
core, having little room to avoid the other hole, resulting 
in large repulsive terms contributing to the anti-binding. 
The electrons, shown in Fig.[5]Cl and C4 are more spread 
out relative to the hole when compared to the binding ex- 
ample of Fig. [5j albeit with a smaller overall volume. 
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FIG. 6. Calculated exciton-exciton interaction energy, 
A xx , against shell thickness for CdS/ZnSe Type-II colloidal 
nanocrystal, with the radial form of the confinement poten- 
tial shown in upper inset. Core diameter of 3.1 nm, with di- 
electric constant of e r =8 (squares with solid line) and for a 
weaker dielectric constant of e,-=5.4 (diamonds and dashed 
line). Lower inset shows A xx against e r , for a fixed shell 
thickness of 2.5 nm. 



V. TYPE II CdS/ZnSe 

There have been several publications focus ed on a simi- 
lar Type-II colloidal core/shell nanost met ure)^2l namely 
CdS/ZnSe. The band gaps and offsets for this particu- 
lar combination result in an opposite potential profile 
compared to a CdTe/CdSe nano-crystal, with the hole 
potential minimum in the shell, and the electron mini- 
mum in the core. A similar transition from binding to 
anti-binding of the biexciton is expected in this Type-II 
structure. 

Particularly noteworthy is the extremely strong mea- 
sured biexciton anti-binding of ^100 me V which led to 
a demonstration of single exciton optical gainP Pertur- 
bation theory results provide some support for these 
values when low values of the dielectric constants are 
employed.^ Again in this case perturbative calculations 
overestimate the anti-binding properties and miss the 
transition from binding to anti-binding. For our model 
system we use a CdS core with a bulk band gap of 
2.485 eV and a ZnSe shell with a bulk band gap of 
2.720 eV, and a conduction band offset of 0.795 eV as 
shown in Fig. [6] (inset) P Electron and hole masses of 
0.2 m e and 0.6 m e respectively are used. We perform all 
simulations at 300 K. 

In Fig. [6] we see a similar transition as in the 
CdSe/CdTe Type II structures, with a transition from 








-30 1 1 ' u 

12 3 

Shell Thickness (nm) 

FIG. 7. Exciton-exciton interaction energy, A xx , plotted 
against increasing shell thickness for ZnSe/CdSe inverted 
Type-I colloidal nanocystal with 3.0 nm (circles with solid 
line) and 3.9 nm (stars with dashed line) core diameters. The 
inset shows the radial form of the confinement potential and 
lines are a guide for the eye. Simulations are performed at 
300 K. 



binding to anti-binding at a shell thickness of around 
0.7 nm, only slightly more than one monolayer shell cov- 
erage. The anti-binding is significantly less than the 
^100 me V measured in Ref. El the uncertainty in di- 
electric constant, masses and core/shell sizes and shapes 
may account for some of this difference, but is unlikely 
to increase the maximum to this level within the shell 
thickness range described here. Fig. [6] also shows the 
effect of a decreased dielectric constant from e r = 8 to 
e r = 5.4, and this is found to give both larger binding 
and anti-binding values. However, a maxima is found 
at around e r = 5.4, with lower dielectric constants than 
this leading to weaker anti-binding or even binding be- 
havior, as shown in Fig. [6] (lower inset). Lower dielectric 
constants lead to an increase in the coulomb interaction 
strength, and in the limit of small dielectric constants the 
four attractive terms can dominate, causing binding. In 
the limit of large dielectric constants, the coulomb inter- 
action strength tends to zero leading to a Axx of zero. 
The maxima lies somewhere in between these two limits, 
here at approximately e r = 5.4. 



VI. INVERTED TYPE-I ZnSe/CdSe 

It has further been suggestecP^ that a similar transi- 
tion in the biexciton binding may be seen in an inverted 
Type-I structure, such as in the ZnSe/CdSe core/shell 
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FIG. 8. Conditional probability densities are shown for a 3.9 nm ZnSe core diameter and 0.25 nm CdSe shell thickness in row 
A. For a 3.9 nm CdTe core diameter and 1.3 nm CdSe shell thickness in row B, and for a 3.9 nm ZnSe core diameter and 3nm 
CdSe shell thickness in row C. The radial form of the confinement potential for each is illustrated. Shown in column 1) is 
g ee , a conditional electron (falling within the blue rectangle) and the resulting electron distribution, column 2) shows g e h, a 
conditional electron and resulting hole distribution, column 3) shows gu e , a conditional hole (falling within the red rectangle) 
and resulting electron distribution and column 4) shows guh, a conditional hole and resulting hole distribution. 



structure, where the electron and hole both have poten- 
tial minima inside the shell. The band offset for the va- 
lence band is small (0.14eV), and for small shell thick- 
ness the hole may delocalize across the entire structure. 
By comparison the conduction band offset is much larger 
(0.86 eV)P Therefore a scenario where the hole is delo- 
calizcd across the structure, and the electron confined to 
the shell, may lead to a quasi Type-II structure.^ As the 
shell thickness is increased we would therefore expect in 
a simple picture to go from binding where both electron 
and hole are localized in the core, to anti-binding as de- 
scribed in the quasi Type-II scenario, and then back to 
binding, where again the electron and hole are localized 
in the shell. 

We use our PI-QMC calculations to go beyond this 
simple picture and include correlations as shown in Fig. [7] 
Compared to the simple discussion above however, we see 
in fact the opposite trend. Due to the correlations we find 
that the biexciton is bound for all shell thicknesses. 

As shown in Fig. [8] sub-panel A1-A4, for a 3.9 nm ZnSe 



core and a thin CdSe shell of 0.25 nm, both the elec- 
trons and holes are confined to the core of the dot, and 
a bound biexciton forms. As the shell thickness is in- 
creased to 1.3 nm (Fig. [8] sub-panel B1-B4), the electron 
and hole (more slowly, due to the lower valence band off- 
set) become more confined to the shell. In Fig.[8]Bl the 
electrons are seen to repel each other, however, not so 
strongly as to localize on the opposite sides of the dot as 
in the previously described Type II anti-binding cases. 
This same feature is evident in the hole-hole repulsion 
shown in B4. At the same time B2 and B3 show excitons 
forming. The bound exciton seen in B2 can be compared 
with the conditional density in B4 where the correlation 
hole matches very closely with the shape of B2, with the 
same feature visible for the electrons in Bl and B3. This 
case can be understood in that there is strong Coulomb 
localization of electrons around the heavier, less well con- 
fined holes. Thus the presence of a hole in the close vicin- 
ity of the electron renders the pair effectively neutral and 
mitigates the inter-electron repulsion. The mechanism 
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FIG. 9. The effect of various parameters on the bind- 
ing to anti-binding regime transition for CdS/CdSe rod/core 
nanocrystal. Red circles with error bars are experimental data 
taken from Sitt et a/.P along with polynomial fit to experi- 
mental data (red dashed line). Electron masses used for each 
dataset are indicated, solid lines indicate a OeV conduction 
band offset. Dashed lines indicated a conduction band offset 
of 0.3 eV. All simulations are with e r =8. Inset shows form 
of confinement potential, black solid line indicates potential 
with OeV conduction band offset, red dashed line shows po- 
tential with 0.3 eV conduction band offset. Lines are a guide 
to the eye. 

for the bound biexciton in this system is that of two in- 
teracting excitons with a weak mutual attraction. As the 
shell thickness increases the increased volume of the dot 
allows like charges to spatially separate more and results 
in a slight increase in binding. 

The binding decreases again for larger shell thick- 
nesses, as the electron and hole become very well confined 
to the shell and like particles are forced closer together, 
increasing the repulsion felt, and reducing the binding 
energy. 

Such complex interplay between particles is handled 
well in the PI-QMC calculation, and as shown, can lead 
to results different to those produced by simple models. 

VII. CdSe/CdS DOT/ROD 

There has also been substantial discussion about col- 
loidal dot /rod structures, in which a colloidal nanocrys- 
tal is embedded inside a nanorod. Sitt et aZP reported 
a similar biexciton binding/anti-binding transition as in 
Type-II nanocrystals. The CdSe/CdS dot/rod struc- 
ture, has been suggested to have either Type-I, or quasi 
Type-II confinement. The debate arises from the cur- 
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FIG. 10. Exciton-exciton interaction energy, A xx , against 
core diameter, for a nanocrystal with a CdS rod and CdSe 
core. Closest fit of PI-QMC simulation (black diamonds) to 
experimental data from Sitt et oZP (red circles, dashed red 
line is polynomial fit to experimental data points) - is with 
OeV conduction band offset and 0.2 m e electron mass, with 
e r =5.785. Inset shows example confinement potential with 
OeV conduction band offset and position of core in rod. 



rent uncertainty regarding the correct conduction band 
offsets. Experimental results however see a strong tran- 
sition from biexciton binding to anti-binding, indicating 
at least some Type-II characteristic behavior. We base 
our model system on that used by Sitt et aZ.P A rod of 
fixed length 40 nm is chosen, and the dot size systemati- 
cally varied. The CdSe core is taken to have a band gap 
of 1.75 eV whilst the CdS rod has a band gap of 2.5 eV, 
and a fixed dielectric constant of 8 is initially used. We 
use 0.3 eV and OeV conduction band offsets to investi- 
gate the effect of these on the exciton-exciton interaction 
energy, A xx . We take the hole mass to be that of the CdS 
core, 0.4 m e , and the electron mass to be either 0.13 m e 
or 0.2 m e . We perform all simulations at 300 K. 

We begin with the core in the middle of the rod. In our 
PI-QMC calculations, Fig.[9j we find that the large offset 
of 0.3 eV with a heavy electron mass, 0.2 m e , confines 
the electron strongly to the core, resulting in an always- 
bound biexciton, resembling Type-I behavior. A lighter 
mass of electron, 0.13 m e , allows it to delocalize further 
outside the dot and along the rod, showing quasi Type-II 
behavior. Similar behavior is seen for a flat offset of OeV 
with the heavier electron mass. A flat offset paired with 
a light electron mass gives a nearly continuously anti- 
bound biexciton, which puts it strongly in the Type-II 
regime. It is clear that the biexciton binding transition 
in the rod is particularly sensitive to the choice of the 



9 



Electron 



Hole 




FIG. 11. Probability densities showing the role of correla- 
tion due to the Coulombic interaction in CdS/CdSe rod/core 
structure, here with a core diameter of 4 nm. The single par- 
ticle electron is highly delocalized along rod, the attractive 
Coulomb potential then strongly localizes the electron to the 
hole, as seen in the exciton and biexciton densities. Black lines 
show the width of the rod (5 nm), with the central 12.8 nm of 
the rod shown out of the total 40 nm length. 



band offset and effective mass parameters. 

In order to match closely the experimental data, we 
use a dielectric constant of the average of CdS and CdSe, 
of 5.785, with a OeV conduction band offset and a heav- 
ier electron mass of 0.2 m P 



As we see in Fig. 10 this 
gives excellent agreement with experimental results. In- 
terestingly, moving the dot nearer to the end of the rod 
was found to have very little impact on the transition. 

We note that the normal perturbative method for cal- 
culating A xx is a particularly poor approximation for 
the dot/rod structures. The single particle electron and 
hole densities for the core/rod structures are shown in 
Fig. [IT] As expected, we see a strongly delocalized elec- 
tron and a localized hole. However, correlation effects in 
these rod structures are striking. Including the Coulomb 
interaction in the path integral calculation, we can see 



the electron density is much more strongly confined to- 
wards the core by its attraction to the hole. This has 
the effect of strongly increasing the binding aspect of the 
biexciton transition. Perturbative methods neglect this 
correlation-enhanced binding; hence, using only pertur- 
bative methods to infer the conduction band offset from 
experimental data can be misleading. 



VIII. CONCLUSIONS 

We have presented for the first time theoretical cal- 
culations for colloidal quantum dots, which systemat- 
ically show the exciton-exciton interaction energy un- 
dergoing a transition from the strongly binding regime 
to the strongly anti-binding regime. Our results illus- 
trate the significance of coulomb correlations in deter- 
mining biexciton binding energies, and show excellent 
agreement with experimental data for CdTe/CdSe Type- 
II nanocrystals over a range of quantum dot sizes. In 
CdS/ZnSe Type-II nanocrystals, we find our results do 
not agree with the large anti-binding values previously 
found, with the effects of higher or lower dielectric con- 
stants unable to account for the discrepancy — raising 
new questions as to the anti-binding mechanism in these 
dots. We also demonstrate that for the inverted Type- 
I nanocrystal system, coulomb correlations can lead to 
a qualitatively different mechanism of binding and anti- 
binding when compared to models which exclude these 
correlations; due to this our results always show bound 
biexcitons. Lastly, we provide insight into novel dot /rod 
structures in which correlations are very strong and per- 
turbative methods are particularly inaccurate. These 
methods will be of importance in assessing the poten- 
tial of particular nanocrystal structures as single exciton 
lasing media. 
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